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Functional integral methods provide a way to define mean-field theories and to sys- 
tematically improve them. For the Hubbard model and similar strong-correlation 
^ . problems, methods based in particular on the Hubbard-Stratonovich transformation 

0^ ' have however been plagued by difficulties to formulate the problem in a spin-rotation 

^•^ . invariant way. Here a formalism circumventing this problem by using a space- and 

time-dependent spin reference axis is discussed. This formulation is then used to 
suggest a possible alternative to Nagaoka ferromagnetism in the strongly correlated 
T^ ' Hubbard model in the vicinity of half-filling. Finally, some aspects of single-particle 

0^ ■ spectra in a simplified model for a short-range ordered antiferromagnet are dis- 

cussed. 

PACS numbers: 71.28. +d, 75.10.Lp, 75.30.Kz 
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^ ■ 1. INTRODUCTION 

o' 

O ' It has been pointed out soon after the discovery of high-temperature super- 

conductivity in the cuprates that electron-electron correlatioas are likely to play 
the crucial role in determining the physics of these compoundsJd This observation is 
based on the fact that compounds with a nominally half-filled band like La2Cu04 

S . are antiferromagnetic insulators rather than metals (as would be expected from 

band theory). This of course is explained by the well-known Mott-Hubbard effect, 
implying important correlation effects. Rather small amounts of doping away from 
half-filling than lead to an unusual metallic state which becomes superconducting 
at high temperatures. _ 

As pointed out by Andersonjd the (strongly correlated) Hubbard model can 
be expected to reproduced many of the essential features: at half-filling the model 
can be mapped onto the antiferromagnetic Heisenberg model and therefore gives 
rise to the antiferromagnetic insulating state, and it is equally clear that doping will 
eventually lead to a conducting state. How exactly this happens is however not clear 
at present: for weak correlation doping initially gives rise to an incommensurate 
antiferromagnetic (or spin density wave) insulator, and only upon further doping 
is a conducting state reachedfl For stronger correlation, the situation may well 
be different. In particular, the spiral state proposed by Shraiman and Siggiatl is 
conducting already at the smallest doping level, even though long-range magnetic 
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order may persist. Moreover, neither the anomalous properties of the normal state 
in the conducting regime nor the origin of superconductivity are currently well- 
understood in the framework of the Hubbard model or of its possible generalizations. 

Whatever the critical value of doping Sc at which magnetic order disappears, 
one can expect that for 5 > S^, short-range antiferromagnetic order with a more or 
less important correlation length ^ still exists. It is then important to understand 
the electronic properties in this regime, in particular in view of understanding the 
unusual properties of the metallic state of the cuprates. Moreover, the ordered state 
at half-filling can be understood in quite some detail (even for strong correlation) 
from Hartree-Fock theory supplemented by spin-wave corrections. At least for 
large (but finite) ^ it therefore seems reasonable to start from a description of the 
ordered state and then, in a second step, to add the effects of the fact that the order 
is only short ranged. Ideally, this might even help to understand the situation of 
very short £,. In fact for large U (and hole doping) a site is either empty and thus 
non-magnetic or singly occupied and thus represents, on some short time scale, 
a magnetic moment. The physics of strong correlation thus is closely related to 
that of more or less short ranged magnetic order. In fact, the dynamics of the 
strongly correlated fermions can be seen as that of fermions moving in a space- and 
time-dependent spin "background" which however is generated by the fermions 
themselves and therefore has to be determined self-consistently. 

In the following section aj^neralization of the standard Hubbard-Stratonovich 
functional integral formalismoEl will be described that explicitly exhibits the local 
spin dynamics and therefore can be expecteito make a formally satisfying treatment 
of the situation discussed above possibleBlJ Some simple applications will then be 
discussed in the subsequent sections. 

2. FORMALISM 

The standard Hubbard Hamiltonian has the form 
H = -tY^iaLa,,, + h.c.) + uY,nr^nri , (1) 

(rr') r 

where aj^ creates a fermion at site r with spin projection cr, i is the nearest neigh- 
bor hopping integral, U is the onsite interaction, and (rr') indicates summation 
over nearest-neighbor bonds, each bond being counted once. Introducing a spinor 
notation via 

the partition function can be written as a functional integral over Grassmann vari- 
ables: 

Z= / X>*e-^»-^-* . (3) 

The free and interaction parts of the action are respectively 



^0-/ dTj^*;(9,-/i)*r-<E(*r*r'+C.C.)l , 
^° [ r (rr) J 



(4) 
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In the following, jthc interaction term will be treated using a Hubbard-Stratonovich 
decompositionfl'la 

For repulsive interactions [U > 0) the appropriate decomposition is 



-eUtp'^l^l^ip^ 
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dAcdAg exp 



U 



(A? 



Here A^ s are real variables, and n = V'tV't + V'tV'i i cr. 
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- leAcU + eAgaz ■ (6) 
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inserts this at each point in space and time and thus obtains a functional integral 
over charge and spin fields Ac_s(''",t), coupled bilinearly to the fermions. A saddle 
point approximation reproduces the Hartree-Fock results, and in particular at half- 
filling one finds an antiferromagnetic (or spin-density wave in another terminology) 
ground state. The unpleasant feature of this way of proceeding is that both Ac and 
As are scalar fields, and one therefore cannot construct easily the effective action 
for the low-energy excitations of the antiferromagnetic state which are spin-waves, 
the existence of which is of course closely related to the vectorial character of the 
order parameter. 

Alternatively, one might use a Hubbard-Stratonovich decomposition using a 
vector auxiliary field. One then however does not even obtain the Hartree-Fock so- 
lution as a saddle point. A number of other, equally unsatisfactory decompositions 
have been discussed in the literature.El In order to obtain a spin-rotation invsiri 
ant effective action for fluctuations around the Hartree-Fock solution, I noticeQ' 
that in writing down the Hamiltonian the choice of the spin quantization axis is 
a priori arbitrary at each lattice site, and in a functional integral formulation can 
also vary in time. I then leave the quantization axis i~lr{T) arbitrary and integrate 
over all possible Jlr ('''), with the appropriate invariant and normalized integration 
measure at each point in space and time. In practice, this is achieved by introduc- 
ing SU{2) rotation matrices Rt{t) at each point of space and time which satisfy 
i?r iT)'^zRt{'^) = ^r(''') • <T- A convenient choice is 



R{Q.) 



cos(i6l) 



(7) 



cos(i0) 
e*'^sin(i6i) 

where 9 and ip are the usual polar angles. I then introduce identities R^ (t)R^{t) — 
1 at the appropriate places in the functional integral and integrate over all config- 
urations J1i.(t). Finally, new spinor variables are introduced via 

$ = i?+* . (8) 

This means that the ^-particles now have their spin along ±r2r(T). The Hubbard 
interaction term is invariant under this transformation, and now the Hubbard- 
Stratonovich transformation can be used in its form ^) without loosing the spin 
excitations which are contained in the functional integral over Jlr(T). This also 
means that a nonzero saddle point value of the spin field As does not necessarily 
imply the existence of magnetic long-range order. For this to occur, the angular 
degrees of freedom have also to be ordered. 
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Given that at least at half-filling Hartree-Fock theory handles local correlations 
rather well even for large U, one can now start by using a saddle point approximation 
for the scalar fields Ac,s- The partition function then becomes 

Z= / Pn P$ X>(5 e"^"^-^"-^^ . (9) 



Here Shf is the action corresponding to the saddle point, and S^ represents the 
coupling between the angular fluctuations and the fermions. Finally Ss represents 
the massive fluctuations of Ac and As around their respective saddle point values, 
and will be neglected in the following. 

For large U an effective action for the spin degrees of freedom and the doped 
carriers can be derived, because in fact arbitrary space-time variations of r2(r,r) 
can be treated. For simplicity, one can then start from a ferromagnetic saddle point 
which is characterized by lower and upper Hubbard band separated by a gap U. 
One then has 

Shf = rdr\j2Kidr-li+ ^Ua,)'^, - t ^ (<I>;<Dr' + cc.) I , (10) 

■Jo [ r ^ (rr) J 

•^0 [ r (rr'> J 

For the case of electron doping, the chemical potential is somewhere in th£ upper 
Hubbard bands, and the lower Hubbard band then can be integrated outJj In this 
way one obtains the effective action for the local spin orientation and particles in 
the upper Hubbard band order by order in t/U. To zeroth order in t/U I find 

(9, - Ai + U)(j), ~ ^(^r(f - COS^r)(l - 0:0r)] 
^ r 

-t^[a(flr,f^r')0:0r'+C.C.]i . (12) 

(rr) J 

Here (j) refers to fermions in the upper Hubbard band, the spin index being omitted, 
and (/?r,^r are the polar angles of Jlr- The coefficients a{i},r,flr') come from the 
expression for the product of two R matrices and are given by 

a(f^r,f^r') = |a|e'''--' =[(l + ^r-^r')/2]^^^exp[ii(f^r,f^r',^)/2] , (13) 

where A(J7i, fl2, fls) is the signed solid angle spanned by the vectors fli,fl2, J^3,c2l 
and z is the unit vector along z. To next order in t/U one recovers the usual t'^ /U 
antiferromagnetic exchange term.cl 

3. A POSSIBLE INSTABILITY OF THE NAGAOKA STATE 

In the absence of particles in the upper Hubbard band, in S^tf only the purely 
imaginary term remains, which is the Berry phase of an isolated spin 1/2, i.e., as 
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expected, the half-filled Hubbard model becomes a collection of independent spins 
for U = oo. Introducing more fermions, two effects occur: (i) the factors {'L^(j>*(j)j.), 
previously introduced by Shankar from semi-phenomenological argumentsjld cancel 
the Berry phase term whenever there is an extra particle on site r, i.e. one is in 
a spin singlet whenever two particles occupy the same site, (ii) the kinetic energy 
term plays a role: in particular, going around an elementary plaquette (1234) the 
lattice curl of the phases Xrr' equals $1234 = [A(r2i, ^2, $^3) + ^(^3, $14, r2i)]/2, 
i.e. there is an effective magnetic field proportional to the solid angle spanned by 
ill, ■■•, ^4- $1234 is the lattice analogue of the familiar winding number density of 
the continuum nonlinear a model. t2l Note that, while the gauge potential in ( [l3| ) 
depends explicitly on z and therefore is not rotational invariant, the physical fluxes 
are. For coplanar configurations $1234 = 0, i.e. the phases can be removed by a 
gauge transformation of the 0's. One then straightforwardly sees that the kinetic 
term is optimized by a_£erromagnetic arrangement of the spins. This is the familiar 
Nagaoka phenomcnon.E3 

One can now ask the question whether non-coplanar configurations of Jlr with 
a nonzero winding number density can lead to an energy lower than the Nagaoka 
state.E^ A configuration giving rise to a spatially constant effective field (i.e. a 
constant lattice curl of Xrr' ) seen by the fermions is shown in fig.|l|. At the mean- 




Fig. 1. A spin arrangement ("texture") giving rise to a constant effective magnetic 
field acting on the fermions. 

field level, one assumes the ilr to be static, and only the first and third terms in 
eq.(|l^) contribute. It is then easy to convince oneself that if this field is such that 
the lowest Landau level is completely filled and all the other Landau level are empty 
the energies of the ferromagnetic and the textured states are nearly identical. More 
precisely, the energy of the textured state is higher than that of the ferromagnetic 
state only due to the fact that now ilr • ^r' < 1, i-e. the global scale of the kinetic 
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energy term is reduced. This leads to an energy difference 

AE « tn^ , (14) 

in favor of the ferromagnet. Here n is the dopant density, i.e. n — 1 corresponds 
to half-filhng. However, one should notice that the fully polarized ferromagnet 
(J7r = z) is an exact eigenstate of the Hamiltonian, whereas the textured state is a 
mean-field trial state, i.e. its mean-field energy is an upper limit to its exact enero- 
In particular, one can calculate the first order quantum fluctuation corrections :lj 
the collective modes have a spectrum Lu{q) ~ 'v\q\-, with a velocity given by 

V = 2tV2Tm . (15) 

This then gives rise to a fluctuation correction to the ground state energy of order 

i.e. of the same order of magnitude but with opposite sign as the correction due 
to the band narrowing, eq.(Q4). To determine which of the two corrections is more 
important, a detailed calculation taking into account short-distance cutoffs would 
be necessary. In the absence of such a calculation one can only point out that a 
textured state is a possible candidate for the destabilization of Nagaoka ferromag- 
netism. 

The existence of a textured state might in fact have interesting consequences: 
(i) there is a large class of spin textures that all give rise to a constant effective 
magnetic field. It is conceivable that this large degeneracy gives rise to sufficiently 
strong quantum fluctuations to destroy long-range magnetic order. The effective 
magnetic field seen by the fermions would however still persist, i.e. the PT symme- 
try breaking inherent in a textured state would survive. This then could leacLto a 
state similar to those proposed in the context of anyon superconductivity.t£tE3 (ii) 
even if long-range order exists the true ground state would be a spin singlet, due to 
the absence of any spontaneous magnetization. This might explaixL the strong vari- 
ation of the total ground state spin observed in finite-size studies .llZl Moreover, with 
twisted boundary conditions, the small-ckjping ground state is in fact found to be 
fully polarized in finite-size calculations .E3 The twist necessary to lead to the polar- 
ized ground state may well be a remnant of the spin texture which (hypothetically) 
is present in the thermodynamic limit. 

4. SINGLE PARTICLE STATES 

It is clearly of importance to understand the single-particle properties in a 
doped and short-range ordered antiferromagnet . Fully understanding this problem, 
everLbased on the simplified action, eq.(n3), and possibly the correction of order 
i/C/Jj is still a formidable unsolved problem. Nevertheless, some insight can be 



gained transforming eqs.(lO) and ( |ll| ) back to the original '^ variables, via eq.(|2). 
One obtains 

^HF + Sn - / dT-{y^^:{dr-li+^Un,-a)^, 
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-tY,{^;^,, + c.c.)\ , (17) 

(rr') J 

i.e. effectively a problem of electrons moving in a space and time-dependent filed 
of constant magnitude but varying orientation. A priori all possible time and space 
dependent configurations Jlr(''") should be summed over. However, it seems likely 
that the dominant configurations will have some degree of short-range antiferro- 
magnetic order. At least in one dimension such a state can be considered as an 
antiferromagnet disordered by defects, and I therefore consider configuration where 
there is antiferromagnetic order and just one defect changing the sign of the order 
parameter. To obtain a calculationally tractable problem, I further consider static 
configurations. In a truely disordered antiferromagnet with some finite correlation 
length, one of course has a finite density of defects (provided these can still be 
defined), but the present highly simplified calculation may still give some helpful 
ideas. In figs.0 and || I show some defect configurations together with their elec- 
tronic spectrum. In particular, in fig.g, local antiferromagnetic order is preserved 
nearly everywhere. One observes that the separation of states into upper and lower 
Hubbard bands remains unless the width of the defect becomes very small (the low- 
est panel in the figure). If one the other hand in the vicinity of the defect there is 
some form of local ferromagnetic order (fig.|3|), one always has states in the middle 
of the band. One can expect that the two rather different types of spectra for lo- 
cal antiferromagnetic or ferromagnetic alignment persist in higher dimensions. The 
two alternative types of spectra (or spectral densities) thus might be observed in 
photoemission spectroscopy. Moreover, configurations like those of fig.^ because of 
their local ferromagnetic order should give rise to an enhanced spin susceptibility. 
In the present calculations the spin structure Jlr was assumed to be static. 
This is probably a reasonable approximation if the electronic excitation energy (as 
measured from the Fermi level) is high compared to a typical frequency of the spin 
dynamics. If on the other hand one is at low excitation energy, i.e. in somewhere in 
the middle of the spectra of figures and if the band is approximately half-filled, 
the spin dynamics will be effectively fast compared to the electronic dynamics. One 
can then speculate that some form of "motional narrowing" will transform the states 
close to the center of the spectra into plane-wave like states, giving rise to Fermi- 
liquid or possibly non-Fermi-liquid like properties. It should also be noted that the 
spectra shown here are symmetric about the horizontal axis, in particular there is 
no "spectral weight transfer" from the upper to the lower Hubbard band. This is 
due to the neglect of the charge interaction contained in the Ss term in eq. IM) . 

5. CONCLUSION 

I have discussed here a way to give a spin-rotation invariant functional integral 
formulation of the strong correlation problem. This formulation makes the role of 
local (in space and time) magnetic order particularly apparent. For half-filling, this 
formulation allows one to recover the mapping to the antiferromagnetic Heisenberg 
model for larsp U, and more generally to obtain the effective spin dynamics for 
arbitrary UBU As discussed here alternatives to Nagaoka ferromagnetism arise at 
strong coupling, and properties of single-particle spectra for different types of local 
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Fig. 2. Single particle states for electrons moving in a magnetic field of constant 
magnitude but varying orientation (see eq.(O)). The spin configurations are shown 
together with the corresponding spectrum and represent a defect in an otherwise 
antiferromagnetically ordered state. Here in all cases approximate antifcrromagnetic 
order is preserved even within the defect. The width of the defect decreases from 
top to bottom. 
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Fig. 3. The same as fig.pl but here there is a strong ferromagnetic component of 
the magnetization at the defects. 
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magnetic order can be discussed. 



generalized so as to include 
One then obtains a formu- 



It should be pointed out that the method can 
the electron-hole symmetry of the Hubbard model 
lation with a fluctuating reference frame in an SU{2) x SU{2) space. This can 
in particular be used to treat particle-particle and particle-hole type instabilities 
on an equal footing in a Ginzburg-Landau like formulation. For example, one 
can include Kanamori-type T-matrix renormalizations into a microscopic formu- 
lation of a Ginzburg-Landau theory of magnetic ordering. Also, the large U limit 
of the present formulation can be relatedQ to the slave-fermion Schwinger boson 
approach. ta 
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